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Abstract - We study charged nanoparticles adsorbed onto surfactant bilayers using small-angle 
scattering of synchrotron radiation. The in-plane interaction of the particles is well described 
by a DLVO component (measured independently in solution) and a repulsive dipolar interaction 
due to the presence of the interface, with an amplitude close to the theoretical prediction. We 
prove that charged nanoparticles at soft interfaces are well described by the classical model of 
Poisson-Boltzmann and van der Waals terms; as a corollary, they do not experience the like- 
charge attraction reported in the literature for some systems of micron-sized spheres at interfaces. 


Introduction. — The distribution of ions close to 
soft interfaces (such as biological membranes) is a long¬ 
standing problem [T], with wide-ranging fundamental and 
practical implications. The variation in ion density as a 
function of the distance to the interface has been thor¬ 
oughly studied, mostly by X-ray scattering techniques 
US], but very little is known about the individual ion- 
ion interaction, although this interaction and the result¬ 
ing ion correlation play an important role in biological 
processes [5]. Such information could also help validate 
and refine certain assumptions made in theoretical and nu¬ 
merical studies of charged systems at the nanoscale, e.g. 
in implicit solvent strategies m- 

One often assumes that the interaction is accurately de¬ 
scribed by the DLVO model [8], which combines the elec¬ 
trostatic and van der Waals interactions, supplemented by 
steric repulsion. While this model has been very success¬ 
ful in explaining the interaction of macroscopic objects, 
its validity at the nanoscale is still in need of direct exper¬ 
imental verification. 

In this context, measuring the correlations between 
highly charged nanoparticles adsorbed at an interface 
could provide edifying results. However, due to their small 
size and often low scattering contrast, very little data ex¬ 
ists for such systems. In the present study, we use small 
(1 nm in diameter) and multivalent inorganic anions ad- 
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sorbed onto bilayers formed by a nonionic surfactant (Fig¬ 
ure [T]) and determine their interaction in the plane of the 
layers. 

Due to the presence of heavy atoms, the structure factor 
of these particles can be measured (even at low concen¬ 
tration and in a disordered state) using the intense X-ray 
beams produced by third-generation synchrotron facilities. 
Using classical concepts from liquid state theory we then 
describe the structure factor in terms of the interaction 
potential and determine the relevant parameters of the 
latter. 



Fig. 1: Schematic representation of charged nanoparticles ad¬ 
sorbed onto a surfactant bilayer (only the top monolayer is 
represented). 
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Materials and methods. 

Preparation. The system under investigation is the 
nonionic lyotropic lamellar phase of the surfactant Brij 30, 
doped with polyoxometalates (POMs): the phospho- 
tungstate anions [PW 12 O 40 ] 3 ~. In previous work pilTO] 
we have shown that the POMs are strongly associated to 
the heads of the surfactant molecules, forming one layer on 
either side of the surfactant membrane. The preparation 
parameters are the overall volume fraction of surfactant <J)g 
and the volume fraction of POMs in the aqueous medium, 
cj>p. Perfectly oriented homeotropic monodomains are pre¬ 
pared inside flat optical capillaries by thermal annealing. 
No salt was added to the samples, so the only available 
counterions are the protons resulting from the (complete) 
dissociation of the phosphotungstic acid. Their water con¬ 
centration varies between 73 and 436 mM, so we can safely 
neglect the possible contaminating ions. 

All chemicals were purchased from Sigma-Aldrich and 
used without any further purification. Brij 30 consists 
mostly of C 12 E0 4 (tetraethylene glycol dodecyl ether), 
along with a smaller amount of homologous C m EO n 
molecules. Its phase diagram is very similar to that of 
pure C 12 E0 4 jllj . At room temperature, it forms a lamel¬ 
lar phase over a wide concentration range, between 25 and 
85 wt%. The surfactant bilayer is 34 A thick [9]. 

The phosphotungstic acid is completely dissociated in 
water, giving [PW 12 O 40 ] 3 ~ polyoxometalates (POMs) and 
H + . The mass fraction of hydration water in the initial 
powder was determined (by drying at 200 °C) as approx¬ 
imately 15%, in good agreement with the 20% value re¬ 
ported in the literature jl2j and was accounted for in the 
preparation. Aqueous POM solutions at a volume frac¬ 
tion (f>p were prepared by dissolving a known quantity of 
powder in ultrapure water, as discussed in detail in a pre¬ 
vious publication [101 . These solutions were then mixed 
with appropriate amounts of surfactant to yield the de¬ 
sired surfactant volume fraction <pg- 

Sample alignment. The POM-doped lamellar phase 
was drawn (by suction with a syringe) into flat optical 
glass capillaries, 0.1 mm thick and 2 mm wide (VitroCom 
Inc., Mountain Lakes, NJ) which were then flame-sealed. 
Using a Mettler FP52 heating stage, the capillaries were 
heated above the lamellar-to-isotropic transition temper¬ 
ature and cooled slowly (0.1 °C/min) down to room tem¬ 
perature. Upon nucleation, the lamellar phase domains 
are aligned with their bilayers parallel to the flat faces 
of the capillary (homeotropic anchoring). These domains 
grow to span the whole width of the capillary |9j. 

SAXS experiments. All data shown in this paper was 
recorded at the ID02 station of the European Synchrotron 
Radiation Facility (ESRF, Grenoble, France), with an in¬ 
cident beam energy of 12.4 keV and at a sample-detector 
distance of 0.885 m. The scattered x-rays were detected 
with a specially developed CCD camera m The flat faces 
of the capillaries were set perpendicular to the x-ray beam 
(normal incidence). The radially averaged intensity I(q) 


is obtained in absolute units; it is subsequently divided by 
the form factor |F/(< 7 )| 2 of a POM (calculated from the 
coordinates of the constituent atoms) and normalized by 
the known density of particles to yield the structure factor 
S(q) over a scattering vector range 0.02 < q < 0.5 A . 

Measurements under oblique incidence were done at 
the bending magnet beamline BM02 (D2AM) [14] of the 
ESRF, at a photon energy of 11 keV and a sample- 
detector distance of 0.270 m. The detector was a CCD 
Peltier-cooled camera (SCX90-1300, from Princeton In¬ 
struments Inc., New Jersey, USA) with a detector size 
of 1340 x 1300 pixels. The data was pre-processed using 
the bm2img software developed at the beamline. Further 
treatment followed closely a method presented in a previ¬ 
ous publication na and is described in the Supplementary 
Material. 

Analysis. For a given interaction potential, the the¬ 
oretical structure factor is calculated iteratively (based 
on the Ornstein-Zernicke equation with the Percus-Yevick 
closure) using the method of Lado [T5][T?] ; details on the 
implementation are given in references [T5JCEE]. With the 
appropriate modifications, this procedure was used for 
modelling the (three-dimensional) interaction of POMs in 
solution as well as for the (two-dimensional) interaction 
of POMs confined within the lamellar phase, as discussed 
below. 


Results. — 

Interaction in aqueous solution. In order to evaluate 
the effect of the interfacial confinement on the interac¬ 
tion between POMs one needs an accurate determination 
of their interaction in solution, in the absence of confine¬ 
ment. We measured the structure factors S(q) for aque¬ 
ous solutions of POMs, with a (j)p between 1 and 6 vol% 
(Fig. [2]). In the framework of the DLVO theory, the in¬ 
teraction outside the hard core consists of an electrostatic 
repulsion and a van der Waals attraction, considered ad¬ 
ditive. The resulting potential is: 


U S oi(r) 


oo for r < 2Rc 

Uei(r) + UvdwM f° r r > 2 R c 


Vei (r) 

KdwM 
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exp(— nr) 
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where Z e s = aZ\, are is the effective charge; Z^ aie = 3 is the 
bare charge of a POM and the renormalization parameter 
a = 0.9 [T9]. l B = 7 A is the Bjerrum length for water, k 
is the reciprocal Debye length (and it depends implicitly 
on the POM concentration). In the van der Waals term 
(written for solid spheres) we defined p = (r/i?o) 2 /2, with 
Ro = 5.2 A the crystallographic radius of the POMs (posi¬ 
tion of the outermost oxygen atom). The only adjustable 
parameters in the model are the effective core radius Rc 
and the Hamaker constant A. The fit yields Rc = 5.53 A 
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and A = 13 kpT, in reasonable agreement with the litera¬ 
ture values for metal oxides: A = 5 — 12.5 [5Dj. 



Fig. 2: Structure factors of POMs in aqueous solution, for 
different volume fractions (pp indicated alongside the curves 
(black), as well as fits using the DLVO potential (JT]) , in red; 
the curves are shifted vertically in steps of 0.2, starting from 
the most concentrated at the bottom. The parameter values 
are indicated in the text. 

No interaction along the bilayer normal. For parti¬ 
cles dispersed into the lamellar phase, one must consider 
the interaction both in-plane and along the layer normal, 
and hence the dependence of the structure factor on both 
q r and q z . Most of the experiments discussed here were 
performed with the incident beam parallel to z (i.e. per¬ 
pendicular to the bilayers and to the flat faces of the capil¬ 
lary), so that the scattering vector is contained within the 
plane of the layers: q = q r (we neglect the curvature of the 
Ewald sphere). The measured S(q) could nevertheless be 
influenced by a possible interaction along z m- In par¬ 
ticular, should such an interaction exist, S(q) would vary 
significantly with the confinement imposed by (ps , which 
imposes the distance between two particle layers facing 
each other across a water space. 

The position of the maximum q m of the structure fac¬ 
tor is sensitive to the interaction. We show this parameter 
in Fig. [Hi as a function of (pp and in Fig. [3jr as a func¬ 
tion of the in-plane density defined as n = 1 fe’ S 1 

with dp = 34 A the bilayer thickness and Vo = 685 A 3 the 
volume occupied by a POM in the crystal. There is no 
simple dependence of q m with (pp (Fig. [Hi), but all the 
values fall onto the same curve when plotted against n 




Fig. 3: Position of the structure factor maximum q m , plotted 
as a function of the volume fraction of particles (j>p (a) and as 
a function of their in-plane density n (b) for different surfac¬ 
tant concentrations (ps (various symbols and colors). In (c) we 
compare structure factors with the same n, but different (pp 

a o — 2 

and (ps ■ The two curves corresponding to n = 0.75 10 A , 
with (ps = 0.4 and 0.5, are shifted upwards by 0.2. 

(Fig. [H>). Moreover, structure factors measured on sam¬ 
ples with the same n but different (pp and (ps share not 
only the maximum position, but also the details of S(q), 
as shown in Fig. [H for two sets of curves with n = 0.75 
and 110 -3 A (the first set is shifted upwards by 0.2). 
We conclude that the relevant parameter for describing 
the interaction is the in-plane density n, and not the bulk 
density (pp\ the lamellar phase can be seen as a stack of 
particle layers (two for each surfactant bilayer) without 
interaction from layer to layer. 

We confirmed this conclusion by measurements where 
the capillary normal is tilted away from the incident beam 
(oblique incidence); this way, the scattering vector q ac¬ 
quires a q z component and one can determine the interac¬ 
tion between particles along z mm- Such an interaction 
is indeed detected for cps = 0.7, but not at (ps up to 0.6, 
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as shown in the Supplementary Material. In the present 
analysis we only include the latter data points. 

Interaction in the plane of the interface. When ad¬ 
sorbed onto the bilayer the particles still experience the 
DLVO interaction, but the presence of the interface in¬ 
duces additional effects. In particular, due to the differ¬ 
ence in ionic strength between the water medium and the 
bilayer, each particle acquires a dipole moment P [22] . 
leading to a power-law repulsion [23J|2Ij : 


V d 


2 P 2 e £b k B T 

4:TTeoe w r 3 k 2 r 3 


( 2 ) 


where the dielectric constants are e w ~ 80 for water and 
eh « 2 for the bilayer j25j. The magnitude of the dipole can 
be estimated as P = f- ytb/^w [23], yielding a prefactor 
B = 2 Z 2 (eb/e w ) « 0.45. We consider that the interaction 
between nanoparticles adsorbed onto the bilayers is given 
by the DLVO term (JTJ) with the same parameter values as 
in the bulk, but corrected for the presence of the image 
charge HSUS and supplemented by the dipolar term ([2]): 


V(r) = 2 — ^— V e \ (r) + Kdw(r) + V d (r) (3) 

The experimental data are shown in Figure 01 for surface 
densities n from 0.17 to 1.9 10 -3 A . We emphasize that 
there is a single fit parameter -the prefactor B in Eq. m 
for all the 16 curves at various concentrations. The best fit 
is obtained for B = 2, somewhat larger than the theoreti¬ 
cally expected value of B = 0.45. To check the significance 
of this discrepancy, we plot in Figure [5] the goodness-of-fit 
function \ 2 , normalized by its minimum value X 2 lin , as a 
function of the fit parameter B. It is difficult to define 
rigorously the uncertainty on the parameter B, since the 
difference between the data and fits in Figure [4] is system¬ 
atic (due in part to imperfect background subtraction), 
rather than statistical. Another source of systematic er¬ 
ror is the approximation inherent in any closure relation 
for the Ornstein-Zernike formula. Redoing the analysis 
with the hypernetted chain closure (widely employed for 
modelling charged systems) yields B = 3, with a similar 
fit quality. The x 2 (B) function being quite flat close to 
the origin for both closures, the experimental result for B 
is not significantly different from the theoretical predic¬ 
tion. On the other hand, the fit quality decreases rapidly 
for negative values of B , ruling out a hypothetical attrac¬ 
tive component. Thus, we do not detect the like-charge 
attraction recently discussed in the case of microspheres 
at the air-water interface [241I27H29] . 

Conclusion. — We infer that the two-dimensional 
structure factor of trivalent anions of nanometer size ad¬ 
sorbed onto nonionic surfactant bilayers is well described 
by the DLVO model and a supplementary dipolar term 
due to the presence of the interface. The water medium 
and the bilayer are treated as continuous dielectrics. This 
approximation is often used in the study of biomolecu- 
lar systems (such as membrane proteins) 30,.3J] but, to 



Fig. 4: Structure factors of POMs adsorbed onto surfactant 
bilayers, for different in-plane number densities n from 0.17 
to 1.9 10” 3 A (indicated alongside each curve), as well as 
best fits using the total potential ©; the curves are shifted 
vertically in steps of 0.2, starting from the most concentrated 
at the bottom. Data corresponding to different confinement 
values <j>s are identified by distinct symbols and colors. 


our knowledge, it has never before been tested experimen¬ 
tally for the interaction between individual ions. The data 
presented here support this assumption within the exper¬ 
imental precision. 
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Fig. 5: Goodness of fit function y 2 normalized by its mini¬ 
mum value Xmin as a function of the parameter B defined in 
Eq. ©, for the Percus-Yevick and hypernetted chain closures 
(solid and dotted line, respectively). The arrow corresponds to 
the theoretical prediction B = 0.45. 
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Supplementary Material — No interaction from 
layer to layer. — 

Experimental. In the main paper we analyzed the 
structure factors S(q) measured in normal incidence and 
concluded that there is no interaction from layer to layer. 

We confirmed this conclusion by measurements where 
the capillary normal is tilted away from the incident beam 
(oblique incidence); this way, the scattering vector q ac¬ 
quires a q z component and one can measure the complete 
structure factor S(q r ,q z ) [TM2T1 . 

Then, the ^-independent component So(q r ) is sub¬ 
tracted. The result should be sinusoidally modulated 
along q z if an interaction is present and essentially zero 
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Fig. 6: Complete structure factor S(q r , q z ) after subtraction of So(q r ) for the sample with rf>s = 0.6 and cj)p = 6%, for different 
tilt angles of the sample (left). Cuts through S along q z at q r = 0.27 A (right). 



Fig. 7: Complete structure factor S(q r ,q z ) after subtraction of So(q r ) for the sample with (f>s = 0.7 and <j>P ~ 6%, for different 
tilt angles of the sample (left). Cuts through S along q z at q r = 0.27 A and sinusoidal fit (right). 
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otherwise m- We present below the results for cj)p = 6%, 
with (f>s = 0.6 and 0.7 (Figures [HI and 0 respectively). 

Within experimental precision, no interaction can be 
detected for (j)g = 0.6, while at cj)s = 0.7 a repulsive in¬ 
teraction (S — S 0 negative for q z =0) is clearly visible. A 
sinusoidal fit yields the distance between the two particle 
layers (facing across a water layer) as 5 = 11.5 A. For the 
analysis presented in the main paper we only include the 
data points without interaction. 

Theoretical. Based on the above data, we can say 
quite confidently that there is no interlayer interaction up 
to (j>s = 0.6; here, we check whether our theoretical model 
for the interaction agrees with this empirical conclusion. 

The comparison is presented in Figure [5] As shown in 
subfigure a), in this configuration the interacting dipoles 
are anti-parallel, so that the corresponding interaction is 
mostly attractive. The distance between the two layers <5 
is measured directly at (j)s = 0.7 (see above) and estimated 
at lower surfactant concentration via the dilution law: it 
increases by almost a factor of two as <fis decreases from 
0.7 to 0.6 (subfigure b). 

Finally, we estimate the interlayer interaction poten¬ 
tial as a function of the in-plane distance V\ (r) for all 
surfactant concentration at the highest particle concen¬ 
tration (f)p = 6% (Figure [5] c). The overall potential is 
positive; there is no hard-core component, since 6 > 2 Rc 
even at the highest For cj>s = 0.4 and 0.5, Vi(r) is 
much lower than kpT. For 4>s = 0.6 it becomes of the 
order of 0.5 but is still very slowly varying. In con¬ 
trast, for ( j>s = 0.7 V \ (r) reaches a maximum of about 
1.5 kpT and varies markedly for r < 2nm. We conclude 
that our model is in good agreement with the experimen¬ 
tal findings, namely no detectable interlayer interaction 
for cj)s = 0.6 (Figure [6| and mild repulsion for <ps = 0.7 
(Figure O. 
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Fig. 8: a) Diagram of the interaction between particles in different layers, across a water space, b) Vertical distance between 
the layers as a function of the surfactant volume fraction (j>s . c) Interlayer interaction potential Vi as a function of the in-plane 
distance r for the highest particle concentration, (ftp = 6%. The top two curves correspond to the experimental situation in 
Figures [6] and [7] 
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